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Abstract 

The discordance between ‘low’ and ‘high’ primordial D/H measurements can be consid¬ 
erably reduced if the analysis of the H+D profiles accounts for the correlated velocity field 
of bulk motion. We have re-estimated the value of D/H at 2 ; = 3.572 towards Q1937”1009 
(‘low’ D/H = (3.3 ±0.3) X 10“^ |p) and at 2; = 0.701 towards Q1718±4807 (‘high’ D/H = 
(2.0 ±0.5) X lO”"^ 1]^) and found that a single D/H value from the range (3.5 —5.2) x 10“^ 
(2cj C.L.) is sufficient to describe both spectra. This result supports homogeneous mod¬ 
els of BBN. The obtained D/H ratio and the measurements of high ^He abundance in 
extra-galactic H H regions Q and ^Li abundance in metal-poor halo stars |§] are consis¬ 
tent, within the errors, with the predictions of standard BBN for 3.6 < ryio < 5.4, which 
corresponds to 0.013 < Db/iioo < 0.020. 


1 Introduction 


The light elements D, ^He, ^He and ^Li are produced in big bang nucleosynthesis (BBN). In the 
standard hot homogeneous BBN model their relative abundances depend on a single parameter 
r] - the ratio between the number of baryons and photons at the epoch of nucleosynthesis (see 
e-g- @)- The value of p is not, however, predicted by the big bang theory. It can only be 
estimated from observations. To £x rj, the measurements of the (D/H)p ratio (the subscript p 
throughout denotes primordial abundance) can be carried out from spectra of distant quasars 
[|^. Such measurements, together with other three light element abundances, provide the 
complete test of the standard BBN model. 

Recent high-redshift D abundances are, however, very discrepant (cf. [Q and |^). It was 
suggested that there may be inhomogeneity in (D/H)p of at least a factor of ten. This implies 
that big bang nucleosynthesis may have occurred inhomogeneously and there may not be 
a unique primordial deuterium abundance as hrst pointed out in p. 

In these proceedings, we show one possibility how to reconcile the observed (D/H)p val¬ 
ues. We demonstrate that the discrepancy between ‘high’ and ‘low’ D/H measurements may 
be overcome in the framework of the line broadening process which accounts for the spatial 
correlations in the intervening velocity held. 



2 The D/H measurements 


We consider two examples of real data which yielded two limiting D/H valnes in the previous 
studies. The hrst one is the absorbing region with observable hydrogen and deuterium Lya 
and Ly/3 lines at z = 3.572 towards Q1937-1009 where a ‘low’ D/H ratio of (3.3 ± 0.3) x 10“® 
was measured in ^j. The second one is the H+D Lya blend with D/H = (2.0 ± 0.5) x 10“^ 
seen at z = 0.701 towards Q1718+4807 [|T^. The latter D/H value is at the upper end of the 
D/H range measured from QSO spectra. 

The foregoing estimations were obtained through the Voigt prohle htting analysis. This 
analysis ignores the large-scale correlations in the velocity held and thus may produce very 
misleading results when applied to the case where the spatial correlations in bulk motion are 
signihcant (see 0, [1], and 0). 

Our method takes into account the correlated structure of the radial velocity held. The 
model supposes a continuous absorbing region of a thickness L exhibiting a mixture of bulk 
motions of diherent types. The gas motion along a given line of sight is described by a huctuat- 
ing (random) velocity held. For the sake of simplicity, we assume a homogeneous (H I) density 
and kinetic temperature Tkin- The velocity held is characterized by its rms amplitude at and 
a correlation length I > 0 (the reader is referred to m for more details). 

To estimate physical parameters and an appropriate velocity held structure along the sight¬ 
line, n(s), we used a reverse Monte Carlo (RMC) technique. Our algorithm requires the def¬ 
inition of a simulation box for the hve physical parameters : N(H I), D/H, Tkin, o't/vth, and 
Ljl (here Vth denotes the thermal width of the hydrogen lines). The velocity, which is the 
continuous random function of the coordinate s, is represented by its sampled values at equally 
spaced intervals As, i.e. by the vector {ui, V 2 , ■ ■ ■, Vk} of the components parallel to the line of 
sight at the spatial points Sj (see H)- 

The obtained results for the z = 3.572 and ^ = 0.701 systems are summarized in Figure 1. 
Upper left panel in this hgure shows conhdence regions in the ‘N(H I)-D/H’ plane for the 
simultaneously htted blue wings of the H-l-D Lya and Ly/3 prohles (Q1937-1009) when the 
other parameters (Tkin = 1-15 x 10"^ K, at/vth = 1.46, and L/l = 3.6) and the corresponding 
conhguration of the velocity held n(s) are hxed. The contours represent 68.3% (innermost), 
95.4%, and 99.7% (outermost) conhdence levels. The hlled circle marks the point of maximum 
likelihood for model (d) from |jn[ where xlnin degree of freedom is 1.18. Lower left panel 
demonstrates the conhdence range for D/H from the second system (Q1718-I-4807). The vertical 
dashed lines restrict the N(H Tjtotai value in accord with |^. The maximum likelihood point 
(hlled circle) corresponds to D/H = 5.4 x 10“®, N(H I) = 1.71x10^'^ cm“^, Tkin = 1-5 x 10^^ K, 

at = 26 km s“^, and L/l =3.5 with Xmin = 0.804. This point is slightly shifted with respect 

,2 


to the best model (e) from |]^ with Xmin — 0.897 (open circle) because in the present study 
a wider simulation box for the physical parameters was adopted. The contours in this panel 
represent 68.3% (inner) and 95.4% conhdence levels. 

These two measurements of D/H in the high-redshift systems reveal no variations of the 
deuterium abundance within the quality of the data sets. We can, therefore, compare SBBN 
predictions with observational abundances which are shown in Figure 1 as well. The theoretical 
SBBN relative abundances of D, ^He, and ^Li (solid curves in the right panels) and their 
uncertainties (dashed curves) as function of rj are taken from |]^ . The abundances of D and ^Li 
are number ratios, whereas Yp is the mass fraction of ^He. In the upper right panel, the rectangle 
height corresponds to the 2cr uncertainty region for D/H from the 2 ; = 3.572 system (which 
provides the highest quality data to date). We do not show here the corresponding bounds 
from the 2 ; = 0.701 system since its H-l-D prohle was obtained with seven times lower signal- 
to-noise ratio resulting in a wider D/H uncertainty range which entirely overlaps the z = 3.572 






















confidence interval for D/H. The rectangle heights in two lower panels give the bonnds from 
recent measnrements of extra-galactic Yp in metal-poor H II regions and ^Li in warm metal- 
poor population II (halo) stars p[. We do not depict an error box representing another ^He 
abundance measurement JT^, since {i) our D/H value is not consistent with Yp = 0.234 ±0.002 
found in 


13 and 


ii 


the estimations of Yp in 


and |13| are not in agreement. The rectangle 


widths give the 2 ct uncertainty regions for rj in the SBBN calculations. The shaded region is 
the window for r] common to all results. It specihes 2 ct conhdence levels in r]. 

Taken at face values, the shaded r^-region 3.6 < rjiQ < 5.4 leads directly to the estimation 
of the present-day baryonic density. Using the present-day number density of photons of the 
cosmic background radiation (T = 2.73 K), we can obtain 0.013 < < 0.020, where 

Hq = lOOh km s“^ Mpc“^. 

Now we can conclude that in the framework of our model we do not detect any variations of 
the D/H ratio. This result supports SBBN. However, more information on the H±D absorption¬ 
line systems is needed to verify this statement unambiguously. 
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Figure 1: Upper left - Confidence regions for the simultaneously fitted H+D Lya and Ly(3 
profiles from the 2 ; = 3.572 system towards Q1937-1009. Lower left - Confidence regions for 
the fitted H+D Lya prohle observed at z = 0.701 towards Q1718+4807. Right - Comparison 
of predicted by SBBN primordial abundances with observational bounds. See text for more 
details. 





































